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Abstract 

An expanding shell is accelerated outward by radiation from the remnant star 
and slowed by ram pressure and accretion as it plows into the interstellar 
medium. We set up the general relativistic equations of motion for such a 
shell. In the non-relativistic limit, the reduce to the standard Ostriker-Gunn 
equations. Furthermore, it is shown that the motion equation can be integrated 
in general, reducing in this way the study of the dynamics of a shell to a set 
of four coupled first order differential equations, solvable up to quadrature. 

Subject headings: Shells, dynamics, explosion. 
PACS number: 97.60L 
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1. INTRODUCTION 



The study of the dynamics of a supernova remnant, a plerion, has been greatly developed 
since the works of Ostriker and Gunn (1971), where they modeled the dynamics of this 
remnant as a thin shell moving in a background which inside has a radiating mass and 
outside is dust. Supposing spherical symmetry, they started with the Newton's second law 
for describing the acceleration of the shell. The force terms were introduced essentially by 
hand. This is the equation of motion: 

where M s is the mass of the shell, M N the mass of the neutron star, R is the radius of the 
shell, P c is the pressure in the cavity, related to the radiation emitted by the star, and Pis 
is the pressure due to the contact with the interstellar medium, defined as Pis = pis (^j) 2 ! 
a friction term known as the "snow plow" effect (Oort 1946). The first term represents the 
gravitational and self-gravitational force, the second a driving force due to the radiation and 
the third a decelerating one due to the interaction with the medium. 

The dynamical study is completed with an equation for the change of the mass of the 
shell (usually due only to the dust that the shell collects while moving in the medium) and 
finally an equation for the change of the shell's internal energy, proportional to the radiation 
emitted by the star. 

This system of equations has been analysed in many ways. Chevalier (1977) and Ostriker 
& Mc Kee (1988) have analysed the kinetic equation of the shell to show the different stages 
of the evolution when each of the terms in turn becomes dominant, while Sato (1988a) and 
Sato & Yamada (1991) have applied this set to the study of the SN1987A in particular. 
Most of the studies use numerical methods to obtain the change in the luminosity or in the 
mass with respect to the change in the radius of the remnant. 



2. SHELLS IN GENERAL RELATIVITY 



On the other hand, the study of the dynamics of a shell separating two backgrounds in 
the context of General Relativity has developed in a powerful and direct formalism since the 
firsts works of Werner Israel (1966), and has been applied to cosmology, mainly to inflation 
(Berezinl987), and to modeling the dynamics of the border between two regions in different 
states, like bubbles in water, or between two given spaces (Sato 1988b). Nevertheless, it has 
not been widely used in the study of the dynamics of the plerion. 

In the context of General Relativity the equation of motion for the shell is obtained 
from Einstein's equations and from the junction conditions for the metric tensor and its first 
derivative. We will work in spherically symmetric spaces, described by 



ds z = -e 2 ^ / dv z - 2be^ dr dv + r 2 dti 2 



(2) 



where / and ip are in general functions of the null coordinate v and of the radius r, and 
b = ±1 depending whether the null coordinate is advanced or retarded. The kinetic equation 
for the shell in this case are the following, with our conventions following MTW (Misner, 
Thorne, & Wheeler 1973), 
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where n a is the 4-vector normal to the shell, u a is the 4-velocity vector , r is the proper 
time, M = 4 7r R 2 a is the proper mass of the shell, a being the surface energy density of the 
shell, and R is the radius. A function in square brackets stands for the difference of that 
function in the region outside, denoted by +, with the region inside the shell, denoted by 
— , and evaluated on the shell. Finally p is the pressure of the shell, and F± is given by 

F± = \ff± + R 2 , (5) 

and dot ' stands for differentiation with respect to the proper time, r. In the general rela- 
tivistic context the dynamical study is completed by the conservation equation, which is an 
equation for the rate of change of the proper mass 

M + pA = AttR 2 [T a pu a n% (6) 

where A = 4-7T R 2 , is the area of the shell. For a derivation of the kinetic equations in some 
particular cases, see for instance (de la Cruz & Israel 1967). 



3. INTEGRATION OF THE MOTION EQUATIONS 



In the present work we want to show that the formalism of shells in General Relativity 
gives not only a motion equation more general than the Newtonian one with tiny correction 
whose measure is "beyond the present technology", but first: includes all the possible terms 
which generate the motion of the shell, no more and no less. Second: even though there are 
terms which in the Newtonian case seem to be highly dissipative, like friction or radiation 
terms, we will show that all these terms can be described within the formalism of General 
Relativity. And Third: we will show that the motion equation, in general, allows a first 
integral, which reduces the problem to a first order dynamical set of equations. We also 



present the Newtonian limits of this integra 
For the line element given by equation ( 
vector, already orthonormalized, at the shel 



of the motion 
]), we have that the 4-velocity and the normal 
are given by 



u a = (v,R,0,0),n a = be^(-R,v,0,0) (7) 
where v is related to R as follows 

The motion equation ([3]) is then rewritten as 

i + M 



where 



/?2 +8 7T P , (9) 



B = Jf + i>, r F+~f, v e^v 2 . (10) 



Now, from the definition of F, equation (EH, we see that 
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so, after some manipulations, the motion equation @ can be expressed as 

The equation for the rate of change of the mass, equation in the spherical case can 
be expressed as 

m = { _h^_ + ^ FR] _ 87l Rp> (13) 
Substituting this last equation in the motion equation (T2), we obtain that 



+ f ) + 



[^^ + lf, v en 2 -^, r FR]=0. (14) 

Finally, using the formula for v, equation (^), can be shown that the expression inside 
the squared brackets is zero!, so 

> + fH0, , 15 ) 

which implies that 

[F] + f = 0, (16) 
the value of the integration constant, C, is zero, as can be shown substituting back [F] + ^ = 



C in the second motion equation (|4j). In this way we have shown that equation ( |T6"D is a 
first integral of the general motion equation of a thin shell with spherical symmetry. 
It proves convenient to introduce a function m = m(u, r) such that 

f± = 1 - (17) 
r 

then, using equation (H), it is easy to show that 

F+ 2 -FJ = - 2 -^, (18) 

where the function m is defined as m = m + — m_, and in some cases can be identified with 
the gravitational massQ. Combining the first integral equation (|T6"D, with equation (0), can 



e 2 

In the cases studied by Israel and de la Cruz or Chase, m± = mi^ — ^r, where mi^ and e\ t 2 
stand for the gravitational masses and electrical charges in both sides of the shell. They defined 
m = m,2 — mi = constant as the gravitational mass of the shell and identified this quantity as the 
conserved total energy of the shell. 
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be obtain an expression for F on either side of the shell in terms of the proper mass and the 
gravitational mass 



m M 

F± = — =F ■ 19 

M 2R K J 

The squared of this last equation , recalling equation ( |T7| ) , results in 

62 /"V m + + m_ M 2 

* " ( M } ~ 1+ i? + 4^' (20) 
which is the motion equation obtain by Lake (1979) from the Lanczos equation. 

4. CONCLUSIONS 

This result of being able to integrate the motion equation in general is quite remarkable 
since it is valid for any matter distribution on either side of the shell, as well as for any 
type of shell. Furthermore, the integration is valid for any way in which the shell interacts 
with the two backgrounds!, the only requirements have been spherical symmetry, the fact 
that the backgrounds satisfy the Einstein's equations, and the conditions of continuity on 
the shell. 

To have a better understanding of the meaning of the first integral of the motion equation, 
it might help to see its form in the Newtonian limit. For this purpose it proves better to 
take the form of the motion equation given by equation (19), considering the cases when 



the masses and the velocity, R, are smaller than the unit and making a Taylor expansion 
we obtain 

1 • M(m_ + f) 

m = M+-MR 2 K -— ^. (21) 

2 R 

The above equation simply states that the sum of the rest energy, the kinetic energy, the 
mutual potential energy and the gravitational self-energy of the shell, respectively in the 
rhs, is not constant, but varies according with the quantity m. 

In this way, the dynamics of a shell moving between any two backgrounds, with spherical 
symmetry reduces to a set of first order equations, namely, the radial motion equation (^0|), 
the balance equation for the proper mass, equation (|13|), the motion equation for the null 
coordinate, equation (H), a state equation for the shell, giving a relation between the pressure 
p and the energy density, a, and for specifying the kind of backgrounds between which the 
shell moves, we have the Einstein's equations, which for this spherical case are (Barrabes & 
Israel 1991) 



m . 



4nr 2 T v r ;m t r = -4 tt r 2 T v v ; i/) >r = 4%rT rr , (22) 



where T^ v is the stress energy tensor, specified by the distribution of matter and energy in 
the chosen background space. 

Despite the above characteristic concerning the dynamics of massive shells, it will be 
worth while to write down conveniently the second order motion equation in order to obtain 
the Newtonian limit, and show how equation ([I]) is recovered. After a direct calculation 
starting from equation (EH), we obtain: 

M(m_ + fF_) rh_ M mF_ M F_ F + 

MR = i -5-2 '- + ^ + —. 23 

R 2 RR R R 
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Note that the first term in the rhs is nothing else than the relativistic generalization of 
the gravitational interaction, and the remaining terms are related with the variation of m_, 
m and M. In fact, we have obtained the relativistic version of the second order equation 
analysed by Ostriker and Gunn (1971) in studying the dynamics of remnants of supernovas 
modelled by thin spherical shells. To show in a clear way the role played by the three last 
terms of the above equation , let us consider that the shell separates two Vaidya spacetimes: 
the interior filled with outgoing radiation and the exterior with ingoing radiation. Thus, 
after taking the Newtonian limit of equation we arrive to the following expression: 



MR= K —- — ^ + 87iRp + 4%R 2 (q_ - q + ) (24) 

R 

where g_ and q + are the energy density of the radiation at each side of the shell and p is 
the pressure of the shell. Therefore, we have recovered, in this simple model, the terms that 
describe the effect of the pressure of radiation inside the cavity and the snow plow effect 
(compare the last term of this equation with the correspondent of the Ostriker's equation ). 
Also, we have an extra term related to the pressure of the shell, since Ostriker has apparently 
considered pressure-less shell. However, a more realistic model is realized in considering the 
exterior spacetime characterized by dust, like in the Friedman- Walker universe, in order to 
taken into account more properly the interstellar medium outside the shell. 

We want to stress the fact that the second order equation , equation (^3|) obtained from 
the first order one, equation (|20"D is completely equivalent to the two second order equations 
obtained from the Einstein's equations, equations ([3], f|) but has a more tractable form and 
the different terms are easier to identify. 

Finally, if instead of the proper time of the shell, r, is used the time of an exterior or 
an interior observer, t> + ,t>_, and then consider the radius of the shell as a function of v, 

R = R(v), so that R = R, v v, after some manipulations can be shown that the motion 
equation for the radius of the shell is given by 

m M 

R, v± = e^A[bA + V Q^^)), (25) 
with rj = ±, for the two possible solutions to the algebraic equation and with 

M R 4 i? 2 y J 

Usually a problem could be posed as follows: given the two backgrounds separated by 
the shell, i. e., two spherically symmetric solutions to the Einstein's equations and a state 
equation for the matter in the shell, determine R, v + , t>_ and M from the motion equations 
which are all first order. In this way we have a set of four first order coupled differential 
equations for 4 unknown functions, so the problem is solvable up to quadrature. 

Nevertheless the analytical integration might prove to be very hard except for the most 
simple cases, so numerical methods are still needed. In a future work (Nunez & Oliveira 
1994) we will present the analysis of this set of equations in various specific backgrounds. 
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